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Abstract 

The smallness of the #13 mixing angle as observed in neutrino oscillation experiments can be 
understood through an approximated g — r exchange symmetry in the neutrino mass matrix. 
Using recent oscillation neutrino data, but assuming no CP violation, we study g — r breaking 
parameter space to establish the conditions under which such a breaking could have a perturbative 
origin. According to the so-obtained conditions, we suggest that a sterile neutrino, matching 
LSND/MiniBooNE neutrino oscillation results, could provide the necessary ingredients to properly 
fix atmospheric and $13 mixing angles to observable values, without exceeding the sterile neutrino 
fraction bound in solar oscillations. In such a scenario, we analyze the general effect of a fourth 
neutrino on the prediction for the effective m ee majorana mass parameter. 
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I. INTRODUCTION 


Neutrino oscillation experiments, using solar, atmospheric, reactor, and accelerator neu¬ 
trinos, provide compelling evidence in favor of nonzero neutrino masses and mixings BE]- 
With the exception of LSND [3], MiniBooNe [l], and a recent reanalysis of the flux in some 
short baseline experiments [5], all existing neutrino oscillation data can be described, and un¬ 
derstood, assuming the mixing of only three flavor (standard) neutrinos. Within this frame¬ 
work, data indicate that two of the three neutrino mass eigenstates, i/ 1>2 , have a squared mass 
difference given by = mj — m\ = A m 2 sol ~ 7.5 x 10~ 5 eD 2 , whereas the third one, z/ 3 , is 
separated from the — u 2 pair by a splitting given by |Am 31 | ~ A m\ TM ~ 2.5 x 10 _3 eU 2 . 
ffowever, the sign in Am 31 = m\ — m 2 , and therefore the neutrino mass hierarchy pattern, 
is still unknown. 

Unlike the quark sector where mixing angles are all small, the measured mixings in 
oscillation experiments are large, except for 0 13 , which has been found to be rather small. 
In the standard parametrization, mixings are given by the PontecorvoMakiNakagawaSakata 
(PMNS) matrix P,E7], 


Upmns — 


C12C13 


S12C13 


sue 


-iScp \ 


-S12C23 — Ci2S23<Sl 3 e 


i&CP 


C12C23 — Si 2 S23"Sl3e 


iScp 


S23C13 


■ K 


( 1 . 1 ) 


V S12S23 - ci 2 C 2 3 si 3 e* 5cp -C12S23 - c 23 si 2 Si 3 e*^ cp c 23 ci 3 


where Cij and Sij stand for cos % and sin 9 t j , respectively, of the mixing angles given as 9 12 , 
0 i3 , and 6 - 23 . Here, 8 cp is the Dirac CP phase, whereas K = Diag(e* /3l / 2 ) e^ 2 / 2 ,1) is a diagonal 
matrix containing two Majorana phases which do not contribute to neutrino oscillations. 
Because of the clear hierarchy in oscillation mass scales, where A rn 2 ATM » A m 2 sol and the 
smallness of $ 13 , it is possible to make a direct identification of above mixings with the 
ones used in a simple two neutrino oscillation analysis of the data. This feature will be 
useful later on for theoretical approximations. Global fits with all three neutrinos indicate 
that [[]], :2] sin 2 d 12 = sin 2 0 Q ~ 0.308 ± 0.017, sin 2 0 23 = sin 2 0 ATM ~ 0 . 437 to;o 2 3 (0.455io;g 33 ), 
and sin 2 0 13 0.0234i t ;o;oo 2 9 ( 0 . 0240 lio;o 82 l)? f° r normal (inverted) hierarchy. S C p, on the 

other hand, has not been determined well so far. 


As in the quark sector, the matrix in Eq. (1.1) actually encodes mixings that are inde¬ 
pendently used to diagonalize both charged and neutral lepton masses. Nevertheless, it is 
always possible to rotate any lepton basis into that where both charged lepton masses and 
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weak interactions are simultaneously diagonal. In such a basis flavor associated to e, //, and 
t, labels became transparent, and, furthermore, the PMNS matrix becomes the one that 
diagonalizes neutrino masses, given in general by the effective operator 

{M l ,) a pP a L(l'/3L) C + h.c ., ( 1 . 2 ) 

such that Upmns — U u ■ K. Therefore, the neutrino mass matrix can be written in terms of 
diagonal (complex) masses, Mdiag = Diaglmie*^ 1 ,77i2e^ 2 , m 3 }, simply as 

M V = U V • M diag ■ U T V . (1.3) 

We will work in such a base hereafter. It is worth noticing that, while the observed 9 13 
is close to zero, although non-null, 9 atm is close to its maximal value, 7 t/ 4. Certainly, 
neither of the central values of these angles is in such critical values; however, it is intriguing 
to observe that, regardless of the hierarchy, it is possible to establish the approximated 
empirical relation 

1/2 - sin 2 9 AT M ~ sin#i 3 /few , (1.4) 

which suggests that the deviation of 9atm from its maximal value, A 9 = 7r/4 — \6 atm\i 
could somehow be correlated to the nonzero value of 0 13 . That would be the case if both 
parameters share the same physical origin. As a matter of fact, in the weak flavor basis we 
have chosen, it is easy to see that null values of A 9 and 0 33 do increase the symmetry in the 
mass neutrino sector, by exhibiting a discrete p—r exchange symmetry [9], Asa consequence, 
observed values of these mixings could be understood as a result of the breaking of p — r 
symmetry. This fact has inspired many theoretical studies in the last years pffT2]. but little 
attention has been paid to exploring models that might provide a physical reason for such a 
breaking. That is the main question we shall address in the present paper by suggesting the 
mixing with a fourth sterile neutrino, that also accounts for LSND/MiniBooNE observed 
oscillations, as the natural source for the violation of p — r symmetry. This idea has been 
explored in Refs. mm , although our general scope in here is quite different. 

The paper is arranged as follows. To clearly establish our sterile neutrino hypothesis, we 
start by revisiting p — r symmetry and parametrizing its breaking. Next, we use experimen¬ 
tal results on neutrino masses and mixings to explore breaking parameter space, assuming 
CP conservation for simplicity, to show that relatively small parameters, and therefore per¬ 
turbative approximations, are well allowed by the data, provided standard neutrino masses 
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are almost degenerate. As we will argue, the order of magnitude of such parameters suggests 
that the naive physical mass scale associated to p — r breaking could straightforwardly be 
identified as the LSND/MiniBooNE scale. Hence, we elaborate a general model for neutrino 
masses and mixings, including a sterile neutrino, and explore the feasibility that the source 
of the breaking came from the sterile neutrino sector, the nonsymmetric couplings of which 
provide for the necessary ingredients to fix all mixings in the model. As we will show, there 
is indeed a non-null region in parameter space where all experimental observables can be 
accommodated within one standard deviation. Furthermore, we calculate the sterile frac¬ 
tion in solar neutrinos predicted by the model and discuss the impact of our sterile neutrino 
model in neutrinoless double beta decay experiments. Finally, we present our conclusions. 


II. NEUTRINO MIXINGS AND /i - r SYMMETRY 


First of all, let us remark that in the theoretical limit of null 0 13 and Qatm — —tt/ 4, with 
only three standard flavor neutrinos, there is not a Dirac CP phase and mixing matrix U v 
becomes the bimaximal mixing form 


Ubm — 


( cos (/?i 2 sin</?i 2 0 

— sin <pi2 cosipi2 —1 

y/2 V2 y/2 

— sin </?i2 cos c/?i2 


( 2 . 1 ) 


Ulli ‘n.z yiz -l 

\ s/2. s/2 s/2 / 

where the only undefined mixing corresponds to 9912 , which eventually, upon small correc¬ 


tions, will become the solar mixing. Using this matrix within Eq. (1.3), one can read out 
the general form of the mass terms, which turn out to be symmetric under the exchange of 
H and r labels. Indeed, by defining the mass matrix elements as m Q a p = (My)^, one obtains 


m a „ = 


mi cos 2 ipi 2 + rn 2 sin 2 ipi 2 ] 


= m °er 


sin 2 ip 


12 


V8 


(m 2 - mi)] 


"V = 




- (mi sin 2 <£12 + m 2 cos 2 tp i2 - m 3 ) ; 
^ (mi sin 2 <£12 + m 2 cos 2 tpi 2 + m 3 ) ; 


( 2 . 2 ) 


where Majorana phases are to be understood. 

Conversely, in the “top-down” approximation, the so-called fi — r symmetry [9] is ex¬ 
pressed as the starting point on mass terms by two general conditions given as m° = m° eT 
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and = mj r , which reduce the number of free mass parameters to 4. Thus, in the limit 
of exact symmetry, one obtains the predictions for mass eigenvalues 

m i = rn° ee — V^m^taiup^, 
m 2 = rn° ee +V2m° eiJ cot (p 12l 

™3 = m - mj r . (2.3) 


where 1 — 2 mixing is given by 

tan2<^ 12 = v^8 


m 


efi 


m 




+ (™°T 


rrvt 


(2.4) 


Besides, null values for d i3 and Ad are predicted. However, as already mentioned, this last 
is not the case from experimental results. Nevertheless, /jl — r can still be assumed as a 
rather approximated symmetry in the neutrino sector, such that understanding the sources 
that contribute to its breaking may enlighten the origin of neutrino mixings. Next, we will 
elaborate on the parametrization for the breaking of /i — t symmetry. 

In general, any generic neutrino mass matrix can always be parametrized in terms of a 
symmetric part plus a correction that explicitly breaks the symmetry, by M u = M^_ T + 5M, 
where M^_ r does posses a /i — r symmetry, whereas 5M is defined by only two nonzero 
elements, 

^0 0 <0 

SM = 0 0 0 , (2.5) 

y 5 0 e ) 

where breaking parameters are, clearly, defined as <5 = m eT — m efi and e = m TT — m Nl . In 
this line of thought, understanding the origin of these parameters is a key to understanding 
di 3 and 9 atm- 


Assuming that these are relatively small parameters, in comparison to m eM and m /J7i 
respectively, observable mixing angles are estimated, in the absence of CP violation, to 
satisfy 


tan2d 0 « 

-/a 

m eil 



V 

1 

+ (m MT - m ee ) _ 

1 


sind 13 ~ 

1 


2 m^S - crrien 

l 

V8 

\_rn% + m^r (■ m M/i - m MT - m ee ) 

sin Ad « 

1 

e (m MT + m ee - m M#1 ) 

i 

3. 

<D 

(M 

+ 


4 

i 

0 ) to 
■ft 

+ 

A 

■c 

1 

■C 

1 

0) 

Cb 
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where m eAt = ( m efl + m eT ) /2 and rli IJ4l = (m MM + m Tr ) /2 . Notice that for small Ad , one 
would have that sin Ad ^ 1/2 — sin 2 0 atm, which jointly to sin0 13 would be given by linear 
relations in terms of e and 5. Of course, the former expressions are first-order calcula¬ 
tions that would provide a good approximation, provided the breaking parameters are small 
enough, ft is remarkable, though, that solar mixing turns out to have a similar expression 
to that obtained in the exact symmetric limit. 

Since we already have quite more precise information about the mixing angles, it seems 
interesting to look at the parameters the other way around, by addressing the theoretical 
question regarding how good p — r is as an approximated symmetry, that is, to obtain 
information about the relative size of the breaking parameters, as a way to search for hints 
of any possible physics lying beneath them. In particular, for instance, knowing to what 
extent 5M could be treated as a perturbation could give a hint toward knowing how far in 
the energy scale the breaking source lies away from the overall active neutrino mass scale. 
This possibility is in itself an interesting one, and our main goal on the following discussion 
will be to explore under which conditions one could achieve a perturbative breaking of p — r, 
meaning acceptable small values for the breaking parameters e and 5. 

Early work has shown that the source of such a breaking cannot come from within the 
Standard Model physics, where the only breaking source is the /i and r mass difference HU- 
As a matter of fact, this charged lepton mass difference is indeed communicated through 
charged weak interactions to the neutrino sector, becoming, upon radiative corrections, a 
source for nonzero 5M. Nevertheless, such a correction turns out to be too small to account 
for observed mixings. Therefore, we are moved to assume that there should be a breaking 
sector out of the Standard Model. 

Without relying on any approximation, one could make a direct reconstruction of the 


mass matrix in Eq. (1.3) and thus of the actual values for <5 and e parameters. To this 
aim, however, one would require knowledge of the mass spectrum, which we do not have 
so far. What we do have, instead, are the two values for mass squared differences involved 
in neutrino oscillations, A m 2 sol and A m? ATM . Thus, one mass parameter in the spectrum, 
which we take as the lightest absolute neutrino mass, aside from the relative sign of mass 
eigenvalues, would remain as free parameters. Notice that by the last we mean to take 
Majorana phases to be either 0 or n so that they provide just a relative sign for the masses. 
Dirac CP phase we will assume hereafter to be zero. In these terms, we rewrite the absolute 
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mass eigenvalues as 


\m 2 \ = \jml + A m 2 sol and \m 3 \ = yjml + \Am 2 ATM \ for NH. (2.7) 

\mi\ = sjml + \Am 2 ATM \ and |m 2 | = \Jml + \Am 2 ATM \ + A m 2 sol for IH. 

Note that, in above, the lightest mass eigenstate, mo, becomes m± for the normal mass 
hierarchy (NH) and m 3 for the inverted mass hierarchy (IH). Next, to proceed with our 
analysis, we define, without approximations, the dimensionless parameters 


5 = 


e = 


= _d_ _ Y,S U ei U Ti-UeiUni) m i 

UeiU^imi 

^(UriUri-U^U^) TYli 

y ^ Ui_iiU 


m, 


e^i 


rn 


( 2 . 8 ) 


W 4 


where the right-hand-sides have been written according to Eq. (1.3). Combined with 


Eq. (2.7), the last expressions give the dimensionless parameters in terms of observed mixing 


angles, oscillation mass scales, and the absolute scale of neutrino masses, mo, as the only 
free parameter. Next, let us perform an approximated analytical analysis of the expressions 


in Eq. (2.8) by considering all four independent combinations of mass signs: (i) m± t 2,3 > 0, 


(ii) mi ;2 < 0 but m 3 > 0, (iii) m lj3 > 0 but m 2 < 0, and (iv) m± < 0 but m 2;3 > 0. Those 
can be written in the suitable form 


5 = 


e = 


y-fsis + y+ 

1 + / Si 3 tan d 2 3 
g cos 2023 ~ si 3 h 

1 + 9 s 23 + s 13^/2 


(2.9) 


where 


V± = 

f = 

9 = 
h = 


C23 ± S23 
C 2 3 

mic \ 2 + crm 2 s ^ 2 — £m 3 

ci 2 si 2 (mi - crm 2 ) 

m i( c i2 5 i3 ~ s i 2 ) + ^m 2 (sl 2 sl 3 - c\ 2 ) + Sw 3 cj 3 
mis \ 2 + am 2 Ci 2 
(mi + <rm 2 ) sin 20 23 sin 26h 2 

m\s\ 2 + crm 2 cf 2 


( 2 . 10 ) 


with the conventions a = +, £ = ± for cases i and ii and cr = —, £ = ± for cases iii and iv, 


respectively. As one can see from these expressions, 5 and e in Eq. (2.9) become zero when 


0\3 = 0 and 023 = —tt/ 4, as expected from exact fi — r symmetry. In the following, let us 
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first examine under which considerations <5 <C 1, and latter on, we will analyze the behavior 
of e in such cases. In the three approaches given by the hierarchies, and using the central 
values for the current mixing parameters, we have 


For NH, rrii <C m 2 ~ \J <C m 3 ~ yjAm 2 ATM , and thus 


/ 


a ^atm I 1 vs 


2 

12 


' , I/I ~ 12.5 , (2.11) 


CIS 12 c 12 \j A m 2 sol ^ S ^ A m\ TM 
which implies |<5| ~ 3.26, discarding NH for any mass sign combinations. 


For IH, mi ~ \J A m\ TM , m 2 ~ ^A m 2 sol + Am\ TM S> m 3 , which gives 

/ 


9 i 9 . crs? 0 

+ usro H-——— 


2 Am 2 ATM 

n r I 1 _ o’ Am L 

S12 C12 1 - (T - o AiiF- 

z LAm ATM 


( 2 . 12 ) 


For cases i and ii, we have |/| rs_/ 10 2 , such that |<5| rs_/ 2, whereas in cases iii and iv, we 
obtain |/| ~ 1 and hence |<5| ~ 0.1. Therefore, cases i and ii again seem to be ruled 
out. 


Finally, for degenerated hierarchy (DH), we get 


/ 


C 12 + VS \ 2 


v S A m 2 A r 

Zj — 77 ——a 


c 2 s 2 
c 12*12 


1 n <7 Am2 sol 

1 “ a 2 m? 


(2.13) 


Cases i and ii give |/| ^ 10 2 , which now implies |h| ^ 15, while for cases iii and iv, one 
gets |/| 1 and |<5| ^ 0.1. Therefore, in DH, cases i and ii are once more disfavored. 


The approximations taken above suggest that only cases iii and iv, in the IH and DH, allow 
for small values of 5. In the IH, a similar analysis, after some algebra, gives e ~ 1, for the 
iii and iv combinations, whereas in the DH, we obtained e ~ 1 for case iii and e ~ 0.4 for 
case iv. Thus, our analysis indicates that the only fairly perturbative case would occur in 
the DH for the signs combination that corresponds to rn\ < 0 and m 2)3 > 0. 

After a complete numerical analysis of the parameter space allowed by data (without 
any approximation), in all four mass sign independent combinations, it was found that, 
while in all possible cases there is always a solution with nonzero values for either of, or 
simultaneously both, the 5 and e parameters, the only case one might consider as fairly 




FIG. 1: One sigma regions for the allowed values of dimensionless /r — r breaking parameters, 
e and 5, as a function of the lightest neutrino mass, m o, for inverted (red line) and normal (blue 
dashed line) hierarchy, in the case where mi < 0 , m 2,3 > 0 . 


perturbative corresponds to mi < 0 and m 2,3 > 0 (regardless of the hierarchy). This is 
consistent with our previous analysis. The allowed one sigma region for both breaking 
parameters in this case is depicted in Fig. [lj Moreover, as one can see from this figure, 
only for almost degenerate neutrinos, where m 0 >0.1 eV, is it possible to actually pick up 
relatively small values for 5 and e to comply with the expectation of a perturbative origin. 
Interestingly enough, none of the breaking parameters is null within such an allowed region. 
We must mention that our results are consistent with those obtained in the general analysis 
made in Ref. [I2j, although our general scope here is quite different. 

Naively, if one takes, for instance, e ~ 5 ~ 0.2 or so, valid in the whole region for mo 


above 0.1 eV for mo, the approximation we used to derive the mixings in Eq. (2.6) would 
be quite well justified, and so would be, to the numerical extent, the geometrical relation 


among mixings given in Eq. (1.4). 


III. fi — r BREAKING FROM A STERILE NEUTRINO 

There are at least two possible approximations one can make to explore the physics be¬ 
yond standard model that is responsible of generating the breaking of /i—r symmetry. Either 
this lies close to the same physics that is responsible for the smallness of neutrino masses, 
in which case one has to probably go for model building to explore concrete possibilities, or 
it is the consequence of the mixing with a sector that does not comply with the symmetry. 
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An example of the latter is the mixing corrections induced through radiative processes and 
due to the explicit violation of the symmetry in the charged lepton masses. As mentioned 
already, this is too small to account for the observed effect in neutrino mixings. Another 
quite straightforward candidate for this would be a sterile neutrino, which by definition 
does not have weak interactions, and thus it has no (e, mu, or tau) lepton flavor. This last 
possibility is much more intriguing, because the not-so-small parameters that are required 
to understand the mixings do suggest that such a new sector cannot be too far away from 
the standard neutrino mass scale. Actually, by assuming that the breaking parameters are 
somehow generated at a given larger scale, m s , and naively taking the perturbations that 
break the symmetry as given in terms of the ratio among the involved scales, which means 
that e , 5 m v /m s , then the mass scale of the sterile neutrino should be just about the eV 
scale, precisely as suggested by LSND/MiniBooNe results. Next, we will analyze in detail 
such a possibility. 

To be specific in our analysis, we assume a single light sterile neutrino and consider its 
most general mass terms, including the mixing with the standard active neutrino sector. 
Notice, however, that we shall be working in a 3 + 1 neutrino mixing scheme in which 
the fourth neutrino (predominantly sterile) is isolated from the block of three active flavor 
neutrinos by the mass gap A m 2 LSND ~ (0.4 — 10) eV 2 . Hence, in the basis (z/ e , z/ M , v T) u s ), 
the mass matrix can be written as 


M 


( M m _ t ~3m s 

l " ct T, m s m s 


(3.1) 


where rn s is the Majorana mass for the sterile neutrino and the vector ~ot T = (a e ,a M ,a T ) 
denotes the active-sterile mixing masses in units of m s . Specific structures of this vector 
could have consequences for model builders, as discussed in Ref. mj. Next, let us assume 
that ag - Cl and m s m a p, where m a p are the elements of the active and explicitly 
symmetric flavor matrix M^_ r . Clearly, if a M = a T , the whole sector would be invariant 
under // — r symmetry, with the known consequences of it for active neutrino mixings. We 
will not assume so, and thus the model will have a single effective parameter for the breaking 
of the symmetry given by the coupling differences Aa = a T — a M . Nevertheless, one would 
find it useful to keep track of the independent a’s along the calculations. 

After decoupling u s , we get, at the lower-order approximation, an effective active flavor 
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matrix, M V) the elements of which are given by the (low-energy) seesaw formula, 


(iWjz)pj — (A/^_ T )p(5 OLpTYl s OLfi . 


(3.2) 


It is clear that M v does not possess in general /i—r symmetry due to the presence of the term 
a p m s aj . It is important to notice that the last can always be separated in a symmetric plus 
a nonsymmetric part under the exchange of p and r indexes. The symmetric part, however, 
will only account for corrections, of second order in a£, to the mass spectrum and the solar 
mixing angle defined by M /1 _ T alone. On the other hand, the nonsymmetric part would be the 
source for the breaking parameters defined in the previous section. It is actually easy to see 
that, without further approximations, one gets 5 = a e Aam s and e = 2a fi Aam s , where, as 
before, = (a p + a T )/ 2. It is straightforward to show that, regardless of hierarchy, our now 
effective dimensionless breaking parameters are second order in ae, and for quasidegenerate 
neutrinos, they can be approximated as 


? \/2m s a e Aa 

o & - , 

m o sin 2p\ 2 

2m s a u Aa 

e « -^- , 3.3 

m 0 cos^ (fi 2 

where the mixing tp 12 is the one involved in the diagonalization of the symmetric sector. 
Within a rough approximation, at lower order, one would have (p 12 ~ Oq. Therefore, to 
get an idea of the order of magnitude of the sterile to active neutrino couplings, one may 
take, for instance, m s ~ 1 eV m 0 ~ 0.2 eV, and e r\_/ r\_/ 0.2 which are consistent with 

the analysis in previous sections, to show that a solution the above formulas is found for 
a e ~ 0.23, a p ~ 0.16, and Aa ~ 0.11. Notice, however, that the effect of ay, even for 
A a = 0, is to incorporate corrections to the mixing in the 1 — 2 sector, and thus, a more 
accurate calculation is likely to modify these naive estimates. 

Notice that, in getting the above results, it seems that three a£ couplings do contribute 
to only two effective breaking parameters, 6 and e. Nevertheless, we would expect that any 
physical solution should at least be around above roughly estimated values for a PtT . To 
address this issue in a more reliable way, one should explicitly confront mass scales and 
mixing angles as obtained by the diagonalization of the complete 3 + 1 neutrino sector 
against measured experimental parameters. To this aim, let us first point out that the mass 


matrix given in Eq. (3.1) contains eight independent parameters (rri£c , a^, and m s ), whereas 
we have knowledge of seven experimentally determined observables, enumerated as follows. 


11 





From weak flavor oscillations one gets two squared mass scales, A m 2 ATM and A m 2 sol , and 
three mixing angles, 9 & , 9 atm , and 6 1 3 . Additionally, from LSND/MiniBooNe results, one 
gets two parameters, taken as a squared mass scale A m 2 LSND ~ A m 2 si ~ rn s and a mixing, 
9 efJ . Therefore, there would be only one free parameter in the analysis, which we take as the 
lightest neutrino mass scale, mo- As already discussed, the consistency of our model with a 
perturbative treatment of the breaking of /i — r symmetry requires m o to be within 0.1 to 
0.4 eV, which corresponds to degenerated hierarchy. This short range for mo will end up 
narrowing the allowed parameter space, as we will show below. 

Next, for our analysis, we will take the intermediate neutrino mass eigenvalues as given 


in terms of the atmospheric and solar scales by Eq. (2.7). Moreover, following the outcome 
of the previous discussion, and considering a perturbative diagonalization of A4, one can 
see that active mass eigenvalues are well approximated (at lower order) by the eigenvalues 
of M^_ r , given in Eq. (2.3), whereas m 4 fa m s . This leave us only with the question of 
constructing a self-consistent system of equations to fit all experimental mixing angles with 
the remaining parameters of the model. By considering the relevant effective oscillations in 
solar, atmospheric, and short baseline experiments, one gets the general formulas 


sin 2 28 & = 4|14,| 2 |Z4 2 | 2 , 
sin 2 2#, 3 = 4|M»3r(|M e i| 2 + |W e2 | 2 ) , 
sin' 20atm — 4|M„ 3 | 2 |«rt| 2 i 
sin 2 2^ ~ 4|M e4 | 2 |M, 4 | 2 , 


(3.4) 

(3.5) 

(3.6) 

(3.7) 


where U ai , for i = 1,2, 3,4, and a — e, /x, r, s, stands for the elements of the general mixing 
matrix which diagonalizes M.. As it is well known, since we are neglecting CP violation, the 
columns of U are given by the properly normalized eigenvectors and A4. Here, we emphasize 


that the left-hand (lhs) sides of Eqs. (3.4 3.6) are known from usual neutrino oscillation 


experimental data, whereas the lhs of Eq. (3.7) comes from considering the allowed regions 
of LSND and Mini-BooNE neutrino data |8j, which we take as 


sin 2 2 6 efl = 0.0023, |Am^| = 0.89eH 2 


(3.8) 


On the other hand, the entries in rhs of Eqs. (3.4 3.7), are given up to 0(a 2 ) by the 
following expressions: 


C e l ~ C12 + Sl 2 " 


ml 


m i2 m i4 


-a_|_a_ — 


c 12 m| 

2 (mu) 2 


a_ 


(3.9) 
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(3.10) 


U 


e2 


Sl2 + C12- 




m 2 i m 24 


-a+a_ 


S12 mi 


2 (m 24 ) 


i a + > 


W, 


e 3 


W e4 


Aa m 2 
V ^2 '^34 Lm 3 i 

m 4 


Ci 2 Si 2 

CK_ + -cy_|_ 

m%2 


QLp 


m 4 i 


«-7= 


1 

V 5 


Aa m| 
2 m 34 


Cl2 Si 2 

-a+- oi— 


m 32 


m 3 1 


+ 


(Aa) 


2 m 2 


4v^ (m 34 ) 2 


U 


/i4 


m 4 

m 4 i 


-a 




U r3 « — 


1 


Aa m 2 r 
2 m 34 


Cl 2 


-a+ 


«12 


-CK_ 


(Aa) 


ml 


4V2 (mu) 2 ' 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


Lm 3 2 m 31 

Here, to simplify, we have introduced the shorthand notation m i} = rrij-mj, for i,j = 1...4, 
a + = a e Si 2 + V / 2ci 2 a /i and a_ = a e Ci 2 — v / 2 si 2 « M , where, as before, c 12 (si 2 ) stands for the 


cosine (sine) function of the free parametric angle, (p 32 , defined by Eq. (2.4). 

As it is easy to see, one can use LSND/MiniBooNe mixing in order to solve for a e in 


terms of a M , using Eq. (3.7). In the quasidegenerate neutrino scenario with sterile mass 


dominance that we are considering, this implies that sin 2d 


e/n 


2| a^a e | 


Numerically, this 


means that (a^ael ~ 0.02 . Similarly, in the same approximation, we obtain for the solar 


mixing 
sin 2 26 P 


o 


sin 2 2<£> 12 


1 - (a 2 + 2a 2 ) + — (cos 2 ip 12 (a 2 e - 2a 2 ) + V^a^ cot 2 <p 12 ) 

Tfl - 


regardless of the hierarchy. It is worth noticing that the last expression does depend on four 
effective parameters, m 0 , m 4 , and a Mjr . In practice, since we are choosing mo in a given 
interval, this relation can be used to formally fix ipi 2 mixing from the equations system, 
leaving us with only two relevant independent parameters: a M and a T couplings. Finally, 
these last parameters can be estimated (at least formally) from the formulas that give 0\ 3 and 
atmospheric mixings, in Eqs. (3.5) and ( |3.6 ), which at the lower order in the a parameters 
are written as 

• 2 on _ ( 4Aam s m 0 \ 2 \ (Am\ TM . 2 2 \ a lx sin 2 (p 12 ] 2 

8111 2013 {~ Aw? a - ' ae * ^~ 2 ± Sm ^ 2 1 ± /x- > ( 3 - 16 ) 


sin 2 26 atm 


(Aa 


ATM 
2 


Am\ TM , _ 2 ) a„sin2»j 12 
4mg ± V 12 I ± ^ 


(3.17) 
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The sign difference in Eq. (3.16) stands for normal and inverted hierarchies, respectively. 
Notice that sin 2 2# 13 comes from corrections at the fourth order in a's, although second 
order in A a, whereas atmospheric mixing gets a second-order correction, as suggested by 


the naive numerical expression in Eq. (1.4). 

Once we have some understanding of the parameter correlations in the determination 


of the four observable mixings given in Eqs. (3.4) to (3.7), we can now proceed with a 
numerical analysis of such a set of equations without further approximations, in order to 
explore and identify the allowed parameter space for a M and a T that gives consistent results 
for current experimental oscillation neutrino data, within one sigma deviations. Our results, 
for mo = 0.2 eV and both the hierarchies, are presented in Fig. [2j where we have scanned 
for appropriated values of a M and a r parameters for each mixing as independent, such that 
the consistent values are found in the overlapping of all regions (shaded area in the given 
plots) that give one sigma values for each standard mixing angle. In these same plots, we 
have also constrained the regions such that —0.4 < e < 0.3 and 0.1 < 5 < 0.6, to insure that 
the whole allowed parameter space be consistent with perturbative approximations. The 
actual effect of including this last condition is to bound the allowed parameter space from 
the bottom and the top, as it can be seen on the given plots. 

By picking up some allowed values for the a parameters, it is easy to reconstruct the 
four-by-four mass matrix to give an explicit numerical example for it. For a typical mass 
matrix obtained by this procedure, we consider 

/ nno^ nm-m nn-Mn moo,) 


M = 


0.0247 0.0110 0.0110 -0.1881 

0.0110 0.0489 -0.0205 0.1315 

0.0110 -0.0205 0.0489 0.0343 

-0.1881 0.1315 0.0343 0.9428 


(3.18) 


for a e = —0.1995, a^ = 0.1395 , and a T = 0.0364 . This matrix leads to the exact active 
neutrino mixings, sin 2 9 & = 0.280 , sin 2 9 atm = 0.379, and sin 2 9 13 = 0 .021, which are 
in good agreement with neutrino oscillation measured parameters, within two, three, and 
two sigma deviations, respectively. Furthermore, we get for the LSND/MiniBooNe mixing 
sin 2 9 e/1 = 0.002, in agreement with the fits of short-baseline neutrino oscillation data [8J. 
On the other hand, the corresponding squared mass differences obtained out of this example 
are A m 2 sol = 7 x 10~ 5 eV 2 and A m\ TM = 2.6 x 10“ 3 eV 2 , whereas we get for the sterile 
mass eigenvalue m s = 0.997 eV, which are also consistent with observations. More accurate 
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FIG. 2: (Color on line) Allowed parameter space for a^ and a T that is consistent with all 

oscillation neutrino data (shaded region), in the 3+1 scenario, within one sigma deviation for 
normal (lhs) and inverted (rhs) hierarchy, respectively. The doted (blue) region corresponds to the 
bounds given by sin 2 20 i 3 5 whereas (red) crossed curves delimit the region for acceptable results 
of sm 2 29 atm- The one sigma sin 2 9 & region is bounded by (black) ex-marked curves. (Yellow) 
squared curves and (gray) star lines bound the regions where e and 5 are small enough for the 
model to be perturbative. 

results could be obtained if higher-order corrections in a parameters are incorporated in 
the reconstruction of the mass matrix. Nevertheless, this numerical matrix serves as a good 
example to illustrate the mechanism we are exploring. 


IV. STERILE IMPACT ON OTHER NEUTRINO OBSERVABLES 


SK and SNO experiments have measured solar electron neutrino flux, $y e , whereas SNO 
has also measured the total solar flux of active neutrinos, Q UetUfl)UT , using neutral current 
interactions. These measurements are in good agreement with the total neutrino flux, d*#, 
predicted by the solar model (see Ref. [2J for further references), which can be used 
to constrain solar neutrino conversion into sterile neutrinos. Thus, assuming that solar 
neutrinos oscillate as u e —> sin av s + cos az/ /t;T , the sterile fraction r\ s = sin 2 a is estimated to 


be [116] 


Vs 




« 0 ± 0.2 . 


(4.1) 
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FIG. 3: (Color online) Allowed ranges of the effective mass \m ee \ as a function of the lightest 

neutrino mass mo in the case of normal (red and light blue) and inverted (green and blue) neutrino 
mass hierarchy. The lower narrowed regions (delimited by the light blue and blue lines) correspond 
to the conventional case of three active neutrinos. Wider upper regions (delimited by red and green 
lines) correspond to our model. 


Theoretically, rj s can be estimated in onr model by using n s = v £* Va —. In terms of mixing 

1 ti/ e —> Ve 

matrix elements, and considering only the contributions at the solar scale, one can write 


Vs 


4UelUe2UslUs2 

i\u a UaY 


(4.2) 


Taking values within the allowed parameter regions for a M and a T , presented in the previous 
section, for ttiq = 0.2 eV, we found that r] s ps (1.2 —1.9) x 10 -2 for normal hierarchy, whereas 
rj s ps (2.7 — 3) x 10~ 2 for inverted hierarchy. Clearly, this results agree with the bounds given 


in Eq. (4.1). 

On the other hand, our sterile neutrino model will also have direct and interesting im¬ 
plications on the effective Majorana. mass term that is involved in neutrinoless double beta 
decay, now written as 


m ee 



(4.3) 


2=1 

The allowed values of \m ee \ in our model can be calculated as a function of the lightest 
neutrino mass. Our results are shown in Fig. [3j where, as before, we have used one sigma 
values for oscillations parameters. Consistently, sterile parameters were considered in the 
regions given as 0.1 < < 0.14, 0.2 < a T < 0.23 for normal hierarchy and 0.12 < < 
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0.14, 0.24 < a T < 0.26 for inverted hierarchy. This ranges are consistent with neutrino 
oscillation data in the region where 0.2 < m 0 < 0.4, which corresponds to the degenerated 
neutrino mass hierarchy. 

As it can be seen, the allowed parameter in our model region is enhanced, compared to 
three-neutrino case, due to the presence of the sterile neutrino. However,However, notice 
that mass hierarchy makes little difference for the allowed parameter space, and thus, it 
would be difficult to be experimentally identified. Nevertheless, other interesting differences 
are at hand. In particular, if forthcoming experiments were to observe a positive signal 
between 0.01-0.4 eV, a degenerated mass spectrum with |m 0 | between 0.2 and 0.4 eV might 
still be possible in the case of four neutrinos. On the other hand, the nonobservation of a 
signal in experiments like GERDA ra would practically rule out this model, which makes 
it falsihable. 

V. CONCLUDING REMARKS 

Neutrino mass models based on fi — r symmetry remain as an interesting possibility 
since they can provide a natural understanding for the almost maximal value of atmospheric 
neutrino mixing, and the smallness of reactor # 13 mixing, using only a couple of generic pa¬ 
rameters that encode the breaking of the symmetry. As the analysis shows, current neutrino 
data are consistent with small values for such parameters, although, in such a scenario, it 
seems to prefer a quasidegenerate active neutrino spectrum. As we have also pointed out, 
the perturbative regime of /i — r symmetry breaking also provides an understanding of the, 
otherwise accidental, relation among atmospheric and 6*13 mixings, which can be expressed 
through the phenomenological numerical formula 1/2 — sin 2 6atm ~ sin# 13 /few. As it turns 
out, from our discussion, at the lower-order approximation, both sides of this equation are 
given as linear expressions in terms of the symmetry breaking parameters. 

On the other hand, the relative smallness of the breaking parameters can, in turn, be 
understood by the mixings of active neutrinos with a sterile neutrino, which by definition 
does not posses a flavor number, and thus neither respects active flavor symmetries. The 
model we have elaborated on in the text incorporates the positive features of (perturbative) 
H — r models, allowing at the same time for a natural explanation of LSND/MiniBooNE 
results. As we have discussed, the model can fix all required parameters using oscillation 
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neutrino observables. The allowed one sigma parameter space turns out to be narrow, 
but we consider it a nice feature of the model, since it allows us to explore its prediction 
without further approximations. In particular, the model is consistent with observed bounds 
on the sterile fraction in solar neutrino flux and predicts distinctive modifications on the 
allowed region for the neutrinoless double beta decay parameter, m ee . From here, EXO limits 
already impose an upper bound for an absolute neutrino mass at about 0.5 eV. Moreover, 
even though the m ee region is wider compared to that of three neutrino scenarios, it predicts 
a lower value for m ee , which would be reachable in forthcoming experiments. As a matter 
of fact, our model could be excluded if no positive signal is found above \m ee \ ~ 0.01 eV. 

Along the analysis, we have not included a Dirac CP-violating phase. Majorana CP 
phases, on the other hand, have been also fixed to 0 or n values, which amount only to 
fixing the relative sing of the mass eigenstates. Nevertheless, it is interesting that our 
analysis shows that the only consistent combination of relative signs that give appropriated 
perturbative solutions comes when m i < 0, whereas other masses are positive. This is an 
intriguing result that could be clarified by an extended exploration of allowed CP phases in 
the model. Such an analysis is out of the scope of the present discussion, but it is part of 
the further work we are already undertaking. 
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